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ABSTRACT 

Line contribution functions are useful diagnostics for studying spectral line formation 
in stellar atmospheres. I derive an expression for the contribution function to the abso¬ 
lute flux depression that emerges from three-dimensional ‘box-in-a-star’ model stellar 
atmospheres. I illustrate the result by comparing the local thermodynamic equilibrium 
(LTE) spectral line formation of the high-excitation permitted OI777 nm lines with 
the non-LTE case. 
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1 INTRODUCTION 

The fundamental parameters of stars such as their effective 
temperature, surface gravity, and chemical composition are 
not observable quantities: rather, they must be inferred us¬ 
ing model stellar atmospheres (Bergemann 2014). Three di¬ 
mensional (3D) hydrodynamic c box-in-a-star’ models (Nord- 
lund 1982) are increasingly being used in this context (Lud¬ 
wig et al. 2009; Magic et al. 2013a; Trampedach et al. 2013). 
These present a huge improvement over classical ID hydro¬ 
static models on account of their ab initio treatment of con¬ 
vective energy transport in the outer envelope that can re¬ 
alistically reproduce the shifting, broadening and strength¬ 
ening of spectral lines by convective velocity fields and at¬ 
mospheric inhomogeneities (Nordlund 1980; Asplund et al. 
1999). Inferred logarithmic abundances can suffer errors as 
large as ±1.0 dex when modelled in ID (Collet et al. 2008). 

Visualizing and understanding spectral line formation 
in three dimensions is non-trivial. Contribution functions 
(de Jager 1952; Gurtovenko et al. 1974) are useful tools 
to that end. They can be interpreted as probability den¬ 
sity functions for line formation in the atmosphere (Staude 
1972; Magain 1986) and are often used to infer the mean 
formation depths of spectral lines. The line intensity contri¬ 
bution function (Magain 1986) represents the contribution 
from different locations of the atmosphere to the depression 
in the normalized intensity. This quantity is commonly used 
to study lines in a solar context (Caffau et al. 2008). Since 
stars are in general not resolved, often more relevant is the 
line flux contribution function, (Albrow & Cottrell 1996), 
which is instead formulated in terms of the depression in 
the absolute flux. 
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Since all parts of the stellar atmosphere contribute to 
its observed flux profile, the line flux contribution function 
is a function of 3D space. Albrow V Cottrell (1996) derive it 
in the context of ID model stellar atmospheres, i.e. assum¬ 
ing plane-parallel symmetry. To apply it directly to a 3D 
model would be to treat the atmosphere as an ensemble of 
ID columns i.e. it would be a 1.5D approximation (Kiselman 
V Nordlund 1995). This is undesirable because the effects 
of horizontal radiative transfer are entirely neglected. An¬ 
other approach is to compute the plane-parallel contribution 
function on a horizontally-averaged, <3D> model. This ap¬ 
proach is still not ideal, because it neglects the effects of the 
atmospheric inhomogeneities which characterize real stellar 
atmospheres. 

In this paper I present in §2 a derivation for the line flux 
contribution function that is valid in three dimensions. To 
illustrate the result, I explore in §3 the formation of the high 
excitation permitted OI 777 nm lines in a 3D hydrodynamic 
STAGGER model atmosphere (Magic et al. 2013a). I present 
a short summary in §4. 

2 THE 3D LINE FLUX CONTRIBUTION 
FUNCTION 

2.1 Concept 

The flux depression at frequency v from a star of radius R 
measured by a distant observer is proportional to the total 
emergent intensity depression, 

Au oc —F j J (/= _ / f ) p dpd(j>, (1) 

in the cylindrical polar coordinate system depicted in 
Figure 1: the polar axis intersects disc-centre and is directed 
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Figure 1. Visual aids to the derivation presented in §2. Left: 
spherical star and an arbitrary reference plane perpendicular to 
the observer’s line of sight; the observer is located to the top 
of the diagram. The spherical polar angle is 6 ; conventionally, 
the notation p, = cos 6 is adopted. The cylindrical polar radius is 
p = R sin 0, where R is the radius of the star, z is the displacement 
from a point in the atmosphere to the arbitrary reference plane. 
Right: spherical star as seen by the observer, who is now located 
out of the page. The azimuthal angle is cj). A ring of constant p is 
shown. 


towards the observer. I v is the specific intensity and II is the 
specific continuum intensity, at some position on the disc, 
in the direction of the observer. The line flux contribution 
function C v must satisfy 

A, = f C„ (r) d 3 r . (2) 

J box 

Crucially, the integration is not performed over vertical 
height as in the plane-parallel derivation of Albrow Cot¬ 
trell (1996), but over the entire 3D volume in which the 
line may form. This is the entire volume of the 3D model 
atmosphere; r thus specifies a position in this box. 

In what follows, equation (1) is manipulated into the 
form of equation (2), and thereby C u is inferred. The con¬ 
tribution function C to the integrated line strength A — 
f Av d v is also found, and satisfies C — f C v dv. 


2.2 Derivation 


Along any given ray, I u and satisfy the respective trans¬ 
port equations, 

(. S „ - /„) , (3) 

(St - It) , (4) 

where z is the path distance, increasing upward towards the 
observer. The linear extinction coefficient a u and the source 
function S u are, in terms of their line and continuum com¬ 
ponents, 


OLv 

Sv 


c , 1 

QLj; H - a v , 

at St + aj sl 

OLv 


(5) 

( 6 ) 


(Hubeny & Mihalas 2014). 

Following Magain (1986), an effective transport equa¬ 
tion for the intensity depression D v = — I v is found by 


subtracting equation (3) from equation (4), 

(Sf - £>„) , (7) 

where the effective source function is 

Sf = Sl) . (8) 

In terms of the optical depth along the ray d r v — — a v d z , 
equation (7) is expressed as 



sf 


(9) 


The formal solution is found by integrating from r u — 0 to 

Tv —> OO, 


Dv 



( 10 ) 


Neglecting proportionality factors, the flux depression is ob¬ 
tained by substituting equation (10) into equation (1), 


Av 



dz p dp dcj ), 


(ii) 


where the integrand is evaluated with the constraint that 
the emergent rays are directed towards the observer. As the 
observer is very far from the star, the emergent rays are 
parallel to each other. Consequently, the last equation is 
written in terms of an infinitesimal volume element, 


Av — l a v St E e Tv d 3 r . (12) 

J star 

The integration in equation (12) is performed over the 
entire volume of the star. 3D box-in-a-star models of stellar 
atmospheres have Cartesian geometry and span a minute 
surface area of the stars they represent (Freytag et al. 2012; 
Magic et al. 2013a). The flux spectrum from the modelled 
star is (approximately) reproduced by shifting the box tan¬ 
gentially across the spherical surface. This is represented by 
two integrations: one over the volume of the box and the 
other over the unit hemisphere. Again neglecting propor¬ 
tionality factors, 

f f a„{r-,Q) Sf (r-,Q) e- T " (r]Q) d 3 rdfi,(13) 

J J box 

where the functional dependence of the integrand has been 
made explicit for clarity. The position vector r specifies a 
position within the box, and the solid angle Q specifies the 
direction of the emergent rays. The infinitesimal solid angle 
satisfies dQ = dfid 4>, where pi — cos 6. After changing the 
order of integration, the contribution function is inferred to 
be, 


Cv (r) = J a v (r; fi) Sf (r; Q) e ~ r ^) d Q . ( 14 ) 

This represents the contribution of a point within the box 
to the observed absolute flux depression in the line, at fre¬ 
quency v. The integrated line strength contribution function 
follows immediately, 

C(r) = J J a v (r; Q) Sf (r;Q) e~ TAr;n) dQ . (15) 


2.3 Rotational broadening 

Line broadening caused by the rigid rotation of the star 
must be included during post-processing. This broadening 
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will affect the monochromatic quantity A v and hence C u . 
Following Dravins &; Nordlund (1990), the broadened spe¬ 
cific intensity is, 

jbroad = g ( ( 16 ) 

where B is a functional which broadens its argument accord¬ 
ing to, 

B [x (u, 0, (j))] = — [ x (v — V sin t sin 0 cos 0 , (ft) d^(17) 

2tt J 

Here u = c^- is the Doppler speed, V is the rotation speed 
of the star in the line forming region, i is the inclination angle 
of the rotation axis with respect to the observer, and the 
integral is over an interval of 271. Retracing the steps above, 
one obtains a rotationally-broadened contribution function, 

C„ (r) = J B [a, (r; Q) Sf (r; Q) e - T -' (r!n) ] dO . (18) 

(In deriving this expression, it is necessary to move B within 
the integral of equation (10). This is valid because the at¬ 
mosphere is assumed to be sufficiently shallow that V does 
not vary across its depth.) 

This integrated line strength A, should not be af¬ 
fected by the rotation of the star (Gray 1992). The adopted 
broadening formalism is consistent with this: integrating 
equation (17) across the line profile, 

J B[x(v,e,4>)] dv = J x (vj 0 , (ft) dv , (19) 

which implies that the contribution function C is not affected 
by the rotation of the star. 


2.4 Mean formation depth 


The interpretation of the contribution function as a proba¬ 
bility density function for line formation (Staude 1972; Ma- 
gain 1986) suggests a formalism for defining the mean for¬ 
mation value of some quantity q with respect to a line, 




Jg(r) C (r) d 3 r 
J C (r) d 3 r 


( 20 ) 


and the variance might then be defined in the usual way as 
E [q 2 ~\ — E [q] 2 . For example, E [q — log 10 r| 00 ] may be used 
to define the mean formation depth, where log 10 t^qq is the 
logarithmic radial optical depth at wavelength A = 500 nm, 
a standard measure of depth in stellar atmospheres. 


2.5 Relationship to the line flux response function 

A related spectral line formation diagnostic is the response 
function: the linear response of the line to a perturbation in 
the atmosphere (Mein 1971; Beckers & Milkey 1975; Caccin 
et al. 1977). The line flux response function 1Z U must satisfy 

5A„ = J Hu (r) 6/3 (r) d 3 r, (21) 

where f3 is an atmospheric parameter (such as temperature). 

Following Magain (1986), the response function is ob¬ 
tained by adapting the above derivation. The effective trans¬ 
port equation equation (7) is perturbed so that D v -A 


D u + S/3 Dft, and the equation for Dft is solved, 

^ = a„ (sf ’ 1 - Dl) , (22) 

where the perturbed effective source function is, 

+ (23) 

The response function is then found by following the pre¬ 
vious derivation, but with D u and replaced by Dft and 
respectively, 

n v (r) = j B [a, (r; fi) Sf ’ 1 (r; fi) e “ T " (r;n) ] dO ,(24) 

and the response function to the integrated line strength is 
7 Z = f IZu dv. 

Response functions can be used to study the sensitivity 
of a spectral line to specific atmospheric variables (Achmad 
et al. 1991). To identify the line forming regions, however, 
contribution functions must be used. 


2.6 Comparison to the plane-parallel line flux 
contribution function 

In the limit of plane-parallel symmetry, the integrand in 
equation (14) loses its dependence on the azimuthal angle (ft: 
a„(r;Q) -4 a v {z;n), 5® ff (r;0) -4 Sf(z;(i), r v (r;Q) -4 
tZ (z; /jl) //x, where z is the geometrical height and rZ is the 
radial optical depth The 3D contribution function C„ thus 
tends to a plane-parallel contribution function C£ p , 

Cr (z) = 2nja v (z; M ) Sf (z; n) d M . (25) 

This expression is the same 1 as that derived by Albrow & 
Cottrell (1996) in the context of ID models, i.e. with the 
implicit assumption of plane-parallel symmetry. 


3 EXAMPLE: 3D NON-LTE SPECTRAL LINE 

FORMATION 

The high excitation permitted OI777 nm lines are known 
to show departures from local thermodynamic equilibrium 
(LTE) (LTE; Sedlmayr 1974; Kiselman & Nordlund 1995; 
Fabbian et al. 2009). It is interesting to consider how the 
lines form within the atmosphere when LTE is imposed, and 
to see what happens once this assumption is relaxed. 

To that end, the contribution function C was im¬ 
plemented into the 3D non-LTE radiative transfer code 
multi3d (Leenaarts & Carlsson 2009). The contribution 
function for the OI 777 nm lines was calculated using a model 
oxygen atom based on those used by Carlsson & Judge 
(1993), Kiselman (1993) and Fabbian et al. (2009). A tem¬ 
poral snapshot of a 3D hydrodynamic model atmosphere 
taken from the STAGGER-grid (Collet et al. 2011; Magic et al. 
2013a) was used. The model was of a typical turn-off star, 
with effective temperature T e ff ~ 6430 K, logarithmic sur¬ 
face gravity (in CGS units) log 10 g — 4, and solar-value 

1 After expressing the contribution function in that paper with 
respect to geometrical height instead of radial optical depth, they 
are the same to a factor of 27t, which arises from those authors 
integrating over spherical polar angle g instead of solid angle Q 
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Figure 2. Material temperature in a vertical slice of a tem¬ 
poral snapshot of a 3D hydrodynamic STAGGER model atmo¬ 
sphere (Magic et al. 2013a). The snapshot has effective temper¬ 
ature T e ff ~ 6430 K, logarithmic surface gravity (in CGS units) 
log 10 9 — A and solar-value abundances. Contours of standard 
logarithmic optical depth log 10 r 5 oo = — 3, — 1, and 1 (from top 
to bottom) are overdrawn. 


abundances (Asplund et al. 2009). The solid angle was sam¬ 
pled using Carlson’s quadrature set A4 (Carlson 1963). 

Figure 2 shows the temperature structure in a vertical 
slice of the snapshot. While the absolute geometrical depth 
and width scales are arbitrary, zero geometrical depth is 
roughly located at the photosphere. Just below this depth is 
the top of the convection zone: hot, light upflows, observed 
as wide granules, turnover to form cool, dense downflows, 
observed as narrow intergranular lanes (Magic et al. 2013a). 
Higher up the atmosphere, reversed granulation patterns can 
be observed: the material above the hot, light granules ex¬ 
pands adiabatically and cools more efficiently, than the ma¬ 
terial above the intergranular lanes - a detailed discussion 
can be found in the appendix of Magic et al. (2013b). 

The LTE and non-LTE contribution functions in this 
snapshot slice are shown in Figure 3. They are both nor¬ 
malized such that the maximum value of the non-LTE con¬ 
tribution function is 1.0. The contribution functions reveal 
that the formation of the lines is qualitatively similar in the 
two cases. There is no contribution at large optical depths. 
This can be attributed the attenuation factor e~ Tu in the ex¬ 
pression for the contribution function, equation (15): deep 
within the atmosphere, photons are more likely to be ab¬ 
sorbed than to penetrate the atmosphere and reach the ob¬ 
server. Line formation is also inefficient in the optically thin 
layers. This is by virtue of the line opacity which appears in 
equation (15): in these layers, a l u ~ 0, so that there is little 
line absorption. Between these two extremes, line formation 
becomes possible once the optical depth becomes small, and 
the factor I c — S\ appearing in the effective source function 
becomes non-zero i.e. once the re-emitted light is no longer 
equal to the absorbed light (Magain 1986; Albrow & Cottrell 
1996). 

Figure 3 shows that imposing LTE inhibits the forma¬ 


tion of the OI777 nm lines. Photon losses in the triplet lines 
themselves (Asplund 2005) drive departures from LTE. The 
line opacity is larger, and the line source function is smaller, 
than their LTE counterparts (Fabbian et al. 2009). This 
leads to a significant strengthening of the lines, correlated 
with the reversed granulation patterns seen in Figure 2. The 
equivalent width ratio is W 3D non - LTE /TY 3D LTE ~ 1 . 5 . 


4 SUMMARY 

Flux profiles observed from stars have contributions from ah 
parts of its atmosphere: thus, the line flux contribution func¬ 
tion is a function of 3D space. In this paper I have shown 
how to derive the contribution function to the absolute flux 
depression that emerges from 3D box-in-a-star model stellar 
atmospheres. The result can be used like other ID contri¬ 
bution functions (Magain 1986; Albrow & Cottrell 1996) to 
help one visualize and understand spectral line formation in 
stellar atmospheres. 
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